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The scattering-based approach for calculating the ballistic conductance of open quantum systems 
is generalized to deal with magnetic transition metals as described by ultrasoft pseudo-potentials. 
As an application we present quantum-mechanical conductance calculations for monatomic Co and 
Ni nanowires with a magnetization reversal. We find that in both Co and Ni nanowires, at the 
Fermi energy, the conductance of d electrons is blocked by a magnetization reversal, while the s 
states (one per spin) are perfectly transmitted, d electrons have a non-vanishing transmission in a 
small energy window below the Fermi level. Here, transmission is larger in Ni than in Co. 

PACS numbers: 73.63.-b, 75.47.Jn, 72.15.-v, 71.15.-m 
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I. INTRODUCTION 

Understanding electron conduction through atomic 
and molecular wires connecting two macroscopic elec- 
trodes is presently a very active research topic. Con- 
ductance measurements for metal nanocontacts display 
flat plateaus and abrupt drops during the elongation at- 
tributed to atomic rearrangements. In noble metals (Cu, 
Ag, Au) and alkali metals (Li, Na, K) the last conduc- 
tance step before breaking, most likely corresponding to a 
monatomic nanocontact or a nanowire, has a value 2Gq 
(G = e 2 /h is the conductance quantum per spin) 
which can be rationalized with the free propagation of 
2 valence s-electrons (one per each spin). For magnetic 
transition metals such as Co and Ni the experimental 
data are less consistent. Oshima et al. 0], who worked 
in vacuum, at variable temperature, and with the pos- 
sibility of an external magnetic field, found conductance 
steps in Ni nanocontacts preferentially near 2Go and 4Go 
at room temperature (RT) and zero field, near 4Go at 770 
K and zero field, and near 3Go (occasionally near Go) at 
RT with a field. Ono et al. Q reported again 2G for Ni 
in zero field, and Go for Ni in a field. Recently, Rodrigues 
et al. Q observed I conductance quantum at RT and zero 
field in a Co atomic chain. Ruitenbeek and collaborators 
H finally obtained low temperature zero field data on 
a variety of magnetic metals, and reported a dominance 
of conductance steps between 2Go and 3Go in Co and 
Ni. Despite the poor agreement between these results 
it is clear that the last conductance step is anomalously 
small with respect to the number of valence electrons de- 
termining the number of conducting channels, since in 
Co and Ni there are many d bands crossing the Fermi 
level. 

Several theoretical methods are available to study the 
transport properties of atomic scale conductors, each 
of them tailored to the underlying electronic structure 



scheme lalalalallialllllialllllllia^ Among these 
are methods based on nonequilibrium Green's functions 
combined with a localized basis set la M, p or with a 
system- independent wavelet basis set |3, [lOJ . A scatter- 
ing approach for conductance calculation has been ap- 
plied by Lang and co-workers ^jjj and Tsukada and co- 
workers ^3] using the jellium model for the electrodes. 
A formulation based on the real-space finite-difference 
approach has been recently proposed and applied to 
calculate conductance of gold and aluminum nanowires 
[l3j | . Recently, the KKR method using the Baranger and 
Stone's formulation of the ballistic transport has been 
presented |L5| . 

In order to deal with magnetic transition metals 
nanowires, we need a formulation capable of handling 
accurately the rather localized 3d electrons. The ap- 
proach to the problem of d electrons based on Vander- 
bilt ultrasoft pseudo-potentials (US-PPs) 17] has gained 
widespread use in recent years. These PPs allow the 
accurate and efficient description of large-scale systems, 
while mantaining all the advantages of a plane waves 
basis set. We already applied US-PPs to several elec- 
tronic structure and nanowire problems [Isj . Recently, 
Choi and Ihm presented a method based on plane 
waves to solve the tip-nanocontact-tip electron scatter- 
ing problem in real atomic contacts and calculate from 
that the ballistic conductance of open quantum systems 
via the Landauer-Biittiker formula. Their approach was 
originally formulated with norm-conserving (NC-PPs) 
Kleinman-Bylander-type pseudo-potentials; but that is 
not very suitable for the localized d electrons. In the 
present paper a generalization of this method to US-PPs 
is presented. As we shall show, it turns out that the 
solution of the scattering problem is not more difficult 
with US-PPs than with NC-PPs, while the calculation 
of the transmission coefficient requires a generalization 
of the expression for the probability current carried by 



2 



a propagating state, which is put forward in this paper. 
The extra workload introduced by the US-PPs is in fact 
negligible. Noting that a method of Thygesen et al. , 
formulated with a wavelet basis set, is also able to deal 
with US-PPs, we will in Sec. IIIA compare our results 
with those of Ref.0 in the specific example of an infinite 
monatomic Al wire with an H atom adsorbed on the side. 

The paper is organized as follows. In Sec. II we de- 
scribe a scattering-based approach for conductance cal- 
culation with ultrasoft pseudo-potentials extending Choi 
and Ihm's method. In Sec. Ill we present some test cal- 
culations and apply our approach to study electron con- 
duction in monatomic magnetic nanowires. A summary 
discussion is given in Sec. IV. 



II. METHOD 

We study electron transport in an open quantum sys- 
tem consisting of a scattering region (0 < z < L) at- 
tached to the left (z < 0) and to the right (z > L) to 
semi-infinite generic electrodes. We assume that the elec- 
trons move ballistically in the self-consistent potential 
with reflection and transmission restricted to the scat- 
tering region. The self-consistent potential is obtained 
performing ground state density functional theory (DFT) 
calculations with a supercell containing the scattering re- 
gion and some portion of left and right electrodes. In 
the xy plane the system is repeated periodically and in 
this plane the scattering states propagating at the en- 
ergy E have the Bloch form and can be classified with a 
k± index. Different do not mix and can be treated 
separately. Furthermore we consider sufficiently large su- 
percells in the xy plane and limit the calculation to the 
two-dimensional T point, k^ = 0. The magnetic prop- 
erties of the system are treated within the local spin- 
density approximation (LSDA). Therefore electrons with 
different spin polarization move independently in differ- 
ent self-consistent potentials and the scattering problem 
is solved separately for the two spin directions. Therefore 
we will from now on solve the problem at some fixed kj^ 
and spin index. Our approach and calculations will ne- 
glect all spin fluctuations, of either quantum or thermal 
origin, assuming a statically spin polarized state of the 
nanocontact. 

We describe transition metals atoms such as Co and 
Ni by ultrasoft pseudopotentials (US-PPs) ^t|- I n that 
scheme the scattering state \1/ at the energy E is a solu- 
tion of single-particle Kohn-Sham equation [2(| (atomic 
units H = 2m =1, e 2 = 2 are used): 



ES\$) = 



-S7 Z + V cS + V NL 



l*>, 



(1) 



|2(j) and Vnl is the nonlocal part of US-PP: 

V N L = J2 Dl mn\PL){P I nl 
I run 



(2) 



constructed using the set of projector functions (3^ as- 
sociated with atom /. The functions (3^ are localized 
within spheres of radius r c centered at the atom /. The 
coefficients D 1 mn characterize the US-PP and depend also 
on the effective potential V e s |20|. The main difference 
with respect to NC-PPs is the appearance in the lhs of 
Eq.Q of the overlap operator: 



S=l + J2lLn\P I m ){P I n\ 



(3) 



as a consequence of the relaxation of norm-conserving 
condition. Here the coefficients q J mn are the inte gral s 
of the augmentation functions defined in Refs. \vA l2fj| . 
Since the energy E is a fixed input parameter one can 
rewrite the Eq.0J as following: 

E*(r) = [-V 2 + Ueff] *(r) + ]T D^(/#|*)/£(r-R,) 

Jmn 

( 4 ) 

with D mn — D mn — Eq I mn , thus recovering a form similar 
to the case of NC-PPs. 

As mentioned above, due to the translational symme- 
try the scattering states have the usual Bloch form in the 
xy plane: 

*(r ± +R ± ,z) = e ik ^V(r x ,z), (5) 
where R_l are the lateral supercell lattice vectors. 

A. Electrode region 

Deep within the two electrodes (z < 0, z > L) the 
scattering state ^ originating from the rightward propa- 
gating Bloch wave i[>j of the left electrode has an asymp- 
totic form: 



* = 



$7 + S 1"ij1pi, Z < 

z > L 



v 3 -r ^ 'ij< 



(6) 



where summation over i G L (i G R) includes the gen- 
eralized Bloch states in the left (right) electrode at the 
energy E which propagate or decay to the left (right). 
The generalized Bloch states (both propagating and de- 
caying) constitute the so-called complex band structure 
of a solid ■ They have the same Bloch form J^J in the 
xy plane while along the z axis, where the translational 
symmetry is absent, they satisfy a condition: 



tpj(r±,z + d) 



ik<! 



(7) 



where V c s is the effective local potential (see Refs.flTl 



where k is in general a complex number and d is the length 
of the unit cell in the z direction of the corresponding 
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bulk crystal. The wave functions ipj with Im/c = and 
Imfc ^ describe propagating and decaying (or growing) 
states, respectively. In the bulk, only the former are of 
course eigenstates. 

Inside the unit cell of the electrode z < z < z + d the 
general solution of the integro-differential Eq.QJ can be 
written as a linear combination |l9j : 

ipj(r) = Y c »J^( r ) + Y c Im,j<t>lm{r)- (8) 

n Im 

Here <p n are linearly independent solutions of the homo- 
geneous equation: 

E<f> n {r) = [-V 2 + V eS (r)] 0„(r), (9) 

and 4>i m is a particular solution of the inhomogeneous 
equation: 

E<j> Im {r) = [-V 2 + 14ff(r)] fam{r) 

+ ^e ik ^(r-r J -Rx), (10) 

where T 1 is the vector-position of the 7-th atom in the 
unit cell. Both <p n and 4>i m are (x, y) periodic as in 
Eq. (jSJ) and the summation over Im in Eq. QHJ involves all 
the projectors in the unit cell. We calculate the functions 
<f> n and 4>i m following the method described in Ref . [Tij . 

The function © is the solution of Eq. (jJJ if the coeffi- 
cients cj mj are given by: 

Clm, j =J2 £)I mn(^j)- (H) 
n 

The allowed values of k for a given energy E can be 
determined by imposing Eq.JJJ) along z to the function 
ipj and to its z-derivative: 

^(r ± ,z + d) = e ik %(r x ,z ), (12) 

^(v ± ,z + d) =e M ^(r ± ,z ). (13) 

Inserting Eq. © into Eqs. (|llll - 113H one can show that 
the last three equations are equivalent to the generalized 
eigenvalue problem: 

AX = e ikd BX. (14) 

where A and B are matrices. We solve this problem 
to obtain in general a complex k and coefficients X — 

|c nj , cimA determining the generalized Bloch state ipj 

at a given energy E and k^. 

B. Scattering region 

Inside the scattering region (0 < z < L) the scattering 
state <J/ can be represented in a form similar to Eq.JSJ): 

*(r) = Y c «^»( r ) + Y cimfamir), (15) 
n Im 



where the functions <p n and 4>im are the solutions of 
Eas. H9I10|> in the region < z < L and the summation 
over Im involves all the projectors within the scattering 
region. The scattering state 'J is completely determined 
by coefficients {c ra , C7 m ,ry ,Uj}. The coefficients cj m are 
now determined by: 

C/m=X^£mO#l*> (16) 

n 

There is one more set of equations on c/ m for the nonlo- 
cal spheres intersecting the boundaries of the scattering 
region and thus shared with the electrodes. One has: 

Cl m ,j + Y r iJ C I™,i = Y Ann(/3nl*> = c Im (17) 
i n 

for spheres intersecting the plane z — and 

Y UjCl m ,i = Y ^In (#1*) = C Jm (18) 
i n 

for spheres intersecting the plane z = L. Ens.l|16 |l -qi8 )l 
together with the usual matching conditions on the wave 
function and its z-derivative on the boundary planes 
z — and z — L give a set of linear algebraic equations. 
Solving this set one obtains the unknown coefficients 
{cniCi m ,fij,tij}. The accuracy of the matching proce- 
dure described above was tested on various systems hav- 
ing no scattering region (e.g., infinite monatomic wires). 
In such systems different propagating modes do not mix 
and must have unit transmission probability (ty = Sij). 
We found this condition to be satisfied with a very high 
accuracy — 5^ | ~ 1CP 6 ). 

C. Conductance calculation 

The ballistic conductance G in the linear regime (in- 
finitely small applied voltage) is related to the total trans- 
mission T at the Fermi energy by the Landauer-Buttiker 
formula G — GqT, where Go = e 2 /h is the conductance 
quantum per spin. The total transmission is given by: 

T = ^|7y| 2 = Tr[T+T], (19) 

ij 

where T is the matrix of normalized transmission am- 
plitudes — y/Ii/Ij ■ tij and Ij is the probability cur- 
rent of the Bloch state ipj in the z direction. Note that 
only rightward propagating states in both left and right 
electrodes should be considered so the matrix T is of di- 
mensions Mr x Ml where Ml and Mr are the number 
of propagating modes in the left and right electrodes, re- 
spectively. The sum over i,j is over the states of both 
spin polarizations present at the chosen energy. In our 
LSDA scheme the two spin channels are decoupled, and 
therefore the Ty is nonvanishing only if i and j have the 
same spin. No spin flip is allowed in our calculation. 
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The eigenvectors of the Hermitian matrix T + T deter- 
mine the coefficients of a unitary transformation from 
the set of Bloch states ipj to the conductance eigenchan- 
nels [2]]]. In the eigenchannel basis the matrix T+T is 
diagonal. Calling Ti the eigenvalues, the conductance is 
a sum of independent contributions from each eigenchan- 
nel: 



G — Go ^ Ti, 



(20) 



where Ti gives the transmission probability for i-th eigen- 
channel. 

In order to calculate the total transmission T one needs 
to know the current Ij carried by the propagating Bloch 
state ipj in the z direction. The ipj are actually pseudo- 
wave functions, and coincide with the all-electron wave 
functions only outside the core regions. Therefore the 
usual expression for the current flowing through a plane 
S perpendicular to the z axis and located at Zq: 



I)' = 21m 



s 



dz 



d 2 r ± 



(21) 



is valid only when the plane z$ does not intersect the 
nonlocal spheres. At a general point it must be modified 
as: 



1 j — 1 i 



21m 



I mn 

f° dz f ^(r-R^'W^n 
J~oo Js 



(22) 



This expression has been derived by requiring the two 
following properties. First, the term added to 1° guaran- 
tees that Ij does not depend on zq once the wave func- 
tion ipj satisfies Eq.(0J and has the Bloch form in the xy 
plane. Second, at all zq where the plane S does not in- 
tersect the nonlocal spheres this additional term vanishes 
so that Ij = Ij. The accuracy of the current calculation 
with Eq. ill'L'l) was tested by checking the unitarity of the 
scattering matrix which is a consequence of current con- 
servation. For each propagating state ipj one must have 
Si + J2i Rij = 1; where TJy and Rij are the normal- 
ized transmission and reflection coefficients, respectively. 
Again, this condition was found to be satisfied with very 
high accuracy for all the systems considered so far. 

We point out that Eg. (|19l) assumes that the current 
operator is diagonal in the basis of Bloch states. That is 
not always true. A mixing might occur between degen- 
erate states with the same k in z direction (for example, 
states derived from atomic p x ,Py or d xz , d yz levels if the 
system has an axial symmetry around z axis). There- 
fore we must orthogonalize Bloch states with respect to 
the current operator before applying Eq. (|19fl for the total 
transmission calculation. The general off-diagonal matrix 



elements Ijk of the current operator are given by: 



Ink - Ij k 

/ dz I /34(r-r J M(r)d 2 rx- (23) 



Iran 



00 
20 



{fi^iY / dz / ^ n (r-T I )]*Mr)d 2 r ± 

) Js 



where 



ip*(r ± ,z ) 
dip](r x ,z] 



dipk(*x,z) 



dz 



(24) 



z=z 



dz 



ipk(r±,z ) 



d 2 r ± . 



The ballistic conductance is calculated in three steps. 
First, we perform the supercell DFT electronic structure 
calculation with the plane-wave code (PWscf ) |22j to ob- 
tain the self-consistent potential V e s and the screened 
coefficients f^„. For spin-polarized systems both V e s 
and T) mn will depend on the spin polarization. Second, 
we calculate the complex band structures of the infinite 
left and right electrodes and orthogonalize the propa- 
gating Bloch states with respect to the current operator 
using Ea. (|23|l . The unit cells of the electrodes are cho- 
sen within the supercell, beyond the scattering region, 
where the potential is bulk-like. Finally, we calculate the 
transmission coefficients ty for each rightward propagat- 
ing state ipj of the left electrode and obtain the total 
transmission T and eigenchannel transmissions T, as de- 
scribed above. This procedure now represents a rigorous 
extension of Choi and Ihm's method to US-PPs. 



III. TEST CASES AND APPLICATION TO 
MONATOMIC MAGNETIC NANOWIRES 

We apply now the above scheme to calculate the bal- 
listic conductance for two test systems: a) a monatomic 
Al wire with an H impurity; b) carbon chains between Al 
electrodes with a finite cross-section. The conductance 
of these systems has been calculated by other methods 
and we can compare our results with the known data. As 
our main application, we study the effect of the magne- 
tization reversal on electron transport in monatomic Co 
and Ni nanowires (a case of ballistic magnetoresistance) . 



A. Monatomic Al wire with an H impurity 

We present first the conductance of a monatomic Al 
wire with an H atom adsorbed on the side and com- 
pare our result with the wavelet basis set calculations 
of Ref . 9] . The supercell for electronic structure calcula- 
tion along with its division onto the bulk-like part (used 
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FIG. 1: Transmission coefficient for monatomic Al wire with 
an H atom adsorbed on the side. The corresponding result of 
Ref. [9] and the transmission coefficient of the perfect infinite 
wire are shown by a dotted and a dashed line, respectively. 



for complex band structure calculations) and the scat- 
tering region are shown in Fig. 1. The LDA calculations 
were performed using the Perdew-Zunger parametriza- 
tion |24j of exchange-correlation energy. The hydrogen 
atom was described by US-PP with the parameters given 
in Ref. 01 . Energy cutoffs of 25 and 250 Ry were used 
for the wave functions and for the charge density, respec- 
tively. The BZ was sampled by three special k-points. 
The spacing between Al atoms was 2.39 A while the wires 
were separated by = 7.5 A. As in Ref.Q the posi- 
tions of the adsorbate and of the two nearest Al atoms 
were optimized while the other Al atoms in the supercell 
were kept at the positions corresponding to the perfect 
infinite wire. The relaxed Al-H distance was found to be 
1.82 A with the Al-H-Al angle of 79.1°. 

In Fig. 1 one can see that the calculated transmission 
curve (solid line) is very similar to that of Ref. |I| (dotted 
line). We recover the main features found in Ref.[j|: a) 
the adsorption of H atom reduces the transmission coef- 
ficient by a factor approximately 0.5 with respect to the 
transmission coefficient of the perfect infinite wire (which 
is merely determined by the number of propagating states 
at a given energy and is shown by a dashed line in the fig- 
ure); b) the drop of the transmission to zero at an energy 
about —2.0 eV. We note a small difference in the positions 
of the edges of transmission curves in two calculations 
(see, for instance, the low maximum at about —1.5 eV) 
which is a consequence of slightly different band struc- 
tures of the bulk- like Al wires. A shift of the energy band 
leads to a corresponding shift of the transmission curve. 



Slight differences of transmission will also be related to 
the optimized atomic positions, presumably slightly dif- 
ferent in our calculations and in those of Ref. [9j . 



B. Carbon chains 

As a second example we consider carbon chains at- 
tached to Al electrodes of finite cross-sections. The elec- 
tron conduction in such systems has been recently calcu- 
lated using methods based on a localized basis set 0, 0] • 
The system consists of a finite length carbon chain placed 
between two Al electrodes of finite cross section oriented 
along Al(100) direction (see Fig. 2). The carbon atoms 
are described by US-PPs with the parameters of Ref. 0] . 
As in the previous case, the energy cutoffs of 25 and 250 
Ry were used for the wave functions and for the charge 
density, respectively, and the BZ was sampled by three 
special k-points. The size of the supercell in the xy plane 
was chosen to be L x = L y = 12.15 A. The electrode unit 
cell is composed of two layers and contains nine atoms. 
The carbon atoms in the chain are separated by 2.5 a.u. 
The carbon chain ends are positioned in the hollow sites 
of the Al(110) electrode sufaces, the electrode-chain dis- 
tance denoted as D. 

We considered seven-atom (with D = 1.0 A) and four- 
atom (with D — 1.9 A) carbon chains previously studied 
with a localized basis set in Ref.[(| and Ref.0, respec- 
tively. Transmission coefficients calculated for these sys- 



scattenng region 




FIG. 2: Transmission coefficient for the seven-atom carbon 
chain between two finite cross-section Al electrodes with the 
electrode-chain separation D = 1.0 A. The dashed line shows 
the corresponding result obtained with a localized basis set 
[6]. 
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FIG. 3: Transmission coefficient for the four-atom carbon 
chain between two finite cross-section Al electrodes as in Fig. 2 
but with an electrode-chain separation D — 1.9 A. The dashed 
line shows the corresponding result obtained with a localized 
basis set [7]. 



tems are presented in Figs. 2 and 3. Our transmission 
curves (solid lines) show qualitativelythe same behaviour 
as the corresponding curves of Refs.jf|,0] (dashed lines). 
Broad maxima in the transmission coefficient are related 
to the energy position of the carbon 7r states while the en- 
ergy windows with vanishing transmission (e.g., energy 
intervals [—1.75,-1.5] and [0.7,1.2]) correspond to gap 
energies where no coupling exists between carbon chain 
and electrode states. We note some small discrepancies 
between the present plane waves calculation and the lo- 
calized basis set calculations, mainly concerning the po- 
sitions of the edges of transmission curves. For example, 
our calculations give zero transmission coefficient in the 
energy window [—1.75, —1.5] to be compared with the en- 
ergy interval [—1.9, —1.7] of Ref.p. As before, the shifts 
of transmission curves appear to be related to slight dif- 
ferences in the band structures of the electrodes obtained 
using two different basis sets. 



C. Monatomic magnetic nanowires 

We now apply the method to calculate the ballistic 
conductance of monatomic Co and Ni nanowires with 
a single spin reversal. We discuss first the electronic 
and geometrical structure of an infinite monatomic Co 
wire without spin reversal |2^] . Similar results for the Ni 
nanowire were reported in Ref. l26l The total energy and 
the magnetic moment of the Co wire as a function of the 
interatomic spacing are shown in the insets of Fig. 4. The 
ground state of the wire is ferromagnetic with an energy 
gain of about 0.41 eV/atom with respect to the nonmag- 



netic wire. The equilibrium interatomic distance for the 
Co wire (about a = 2.11 A) is similar to that of the Ni 
wire [2^1 while the magnetic moment is more than twice 
(/j = 2.16 /xe/atom compared to 1.11 /XB/atom for the Ni 
wire). Similar to the Ni wire |2rj| the magnetic moment 
of the Co wire vanishes and the wire turns nonmagnetic 
in the unphysical region where the interatomic distance 
is taken less than approximately 1.6 A (1.9 A for Ni). 

For a perfect infinite monatomic wire the number of 
conducting channels is just the number N of bands cross- 
ing the Fermi level. Since in the absence of tip-wire junc- 
tion each channel has a unit transmission probability the 
ballistic conductance as given by the Landauer-Buttiker 
formula is G = GqN. In Fig. 4 we show the complex 
band structure of a monatomic Co wire obtained at the 
equilibrium interatomic spacing - an ingredient that will 
be needed for further conductance calculation. Com- 
plex bands for the majority (spin up) and minority (spin 
down) spin polarizations are displayed in Figs. 4a and 4b 
respectively. Qualitatively, the complex bands of the Co 
wire look quite similar to those for a Ni wire [13] though 
the exchange splitting is twice as large (approximately 
2 eV versus 1 eV for the Ni wire) which is directly re- 
lated to a stronger Hund's rule coupling and a higher 
magnetic moment per atom of the Co wire. As for the 
number of conductance channels, from the figure one can 
count one spin up and six spin down (real) bands cross- 
ing the Fermi level, giving altogether seven propagating 
channels. The same result was previously found for the 
Ni wire, jl^, I22} Not surprisingly, this overall number 
of channels is much larger than any of the observed ex- 
perimental conductance steps, of both nickel and cobalt 
nanocontacts 0, IS El HI • 

We consider now a monatomic nanowire with posi- 
tive magnetization as our left electrode, and another 
with negative magnetization as our right electrode, the 
nanocontact thus consisting of a magnetization reversal. 
We obtain the optimal shape of magnetization reversal 
by LSDA calculations, minimizing the total energy of a 
supercell containing 12 atoms. The detailed shape of the 
optimal magnetization reversal was obtained by allow- 
ing magnetization to vary freely in magnitude and sign, 
while constraining the overall total magnetic moment to 
be zero (for computational details see Ref. Er^. The re- 
sulting spin reversal turned out to be quite sharp and the 
ground state had, roughly speaking, 6 atoms with posi- 
tive and 6 atoms with negative magnetization (see Fig. 5) 
While this way to build a collinear magnetization rever- 
sal would be of course incorrect for a bulk Bloch wall, 
it may not be an unreasonable description of magneti- 
zation reversal inside a very short nanocontact. Bruno 
[28| showed that the width of a nanoconstrained mag- 
netic domain wall is essentially determined by the size 
of the constriction and can thus be much smaller than 
the (noncollinear) Bloch wall in 3D crystals or Neel walls 
in a film. Essentially, in such an atomically constrained 
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FIG. 4: Complex band structure of monatomic Co nanowire for the majority (spin up) and minority (spin down) states. Real 
bands, complex bands with Re k — 0, and it /a are plotted in the middle, left, and right panels, respectively. Each band is 
labeled by its main atomic character. The total energy and the magnetic moment per atom as a function of interatomic spacing 
are shown on the left and right insets, respectively. 



system, the domain wall spin rotation must take place in 
such small distance to justify a spin collinear approxima- 
tion. 

The self-consistent potential for the nanowire with a 
spin reversal is subsequently used to solve the scattering 
problem. The scattering region and the unit cell of the 
left and right side of the wire are shown in Fig. 6. For 
symmetry reasons the total transmission is the same for 
electrons of both spin directions (this conclusion has also 
been checked numerically). Therefore we choose to con- 
sider only electrons with spin down polarization incident 
from the left on the spin reversal region. 



i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — r 




5 10 15 20 25 



z (A) 

FIG. 5: Planar average of the magnetization along the Co 
nanowire. 



In Fig. 6a we present the total transmission (for one 
spin channel) and the contributions from each band as 
a function of energy for the Co wire with a spin rever- 
sal. The propagating state in the left half of wire can 
propagate through the spin reversal and contribute to 
the total transmission only if it can be matched with an 
appropriate propagating state to the right. For exam- 
ple, at the Fermi energy there are 6 propagating states 
in the left side. However, only s — d z 2 and d z 2 — s states 
have nonzero transmission since there is only one d z 2 — s 
propagating state in the right side. The d states in par- 
ticular can be matched only with decaying states in the 
right side and are therefore completely reflected due to 
the spin reversal. At energy about E = Ep — 0.4 eV 
the two d xz ,yz states start contributing to the transmis- 
sion because there appear two states of the same sym- 
metry in the right side of the wire. On the contrary, the 
d X y,x 2 -y 2 states cannot get through the spin reversal at 
any energy and never contribute to conductance. Fig. 6b 
shows the eigenchannel transmissions Tj. The s — d z i 
and d z 2 — s states are combined to form a single per- 
fectly propagating conductance eigenchannel (essentially 
s like). Another channel (mostly d z 2), orthogonal to the 
first one, appears at energy about 0.6 eV below the Fermi 
energy. It is quite narrow and has a smaller transmis- 
sion probability (about 0.6 at the maximum). The two 
d X z,yz states form two eigenchannels of the same symme- 
try. They contribute to the total transmission between 
0.4 eV and 1.8 eV below the Fermi energy and are par- 
tially reflected (the transmission probability is about 0.8 
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FIG. 6: Transmission coefficients for spin down electron incident from the left on a magnetization reversal in a monatomic 
Co nanowire: a) total transmission (solid curve) and contributions from each band; b) eigenchannel transmissions T,:. The 
corresponding data for a monatomic Ni nanowire are shown on the insets. The bands and conductance channels are labeled by 
their main atomic character. The Fermi energy is chosen to be zero (dashed vertical lines). 



at the maximum). 

In the left and right insets of Fig. 6 we show the to- 
tal transmission (for one spin channel) and eigenchannel 
transmissions, respectively, similarly calculated for a Ni 
nanowire with a spin reversal. One can see that again 
there is a single conductance channel perfectly propa- 
gating at all the energies studied contributing one con- 
ductance quantum (per spin) to the total conductance 
at the Fermi energy. However, the contribution from 
other d channels at lower energies is now much larger 
(at the maxima these channels have almost unit trans- 
mission probabilities). This seems reasonable since the 
Co atoms in the wire have a larger exchange splitting, 
and are therefore more reflecting, compared to the Ni 
atoms (as discussed above). 

Adding the transmission of the two spin channels at 
the Fermi energy we get the overall value of Go -I- Go for 
the conductance. This is to be compared with Go + 6Go 
of the ferromagnetic wire. Such a large difference be- 
tween conductances of uniform ferromagnetic wire and 
the same ferromagnetic wire with a spin reversal is an ex- 
ample of magnetoresistance, analogous to that discussed 
earlier. |2^| Note that the minimal conductance value 2Go 
obtained in presence of spin reversal is still larger than 



that of "fractional peak" at approximately Go reported 
for Co and Ni nanowires 0. We should stress, how- 
ever, that the exact equivalence of the s channels for 
the two spin polarizations, leading to the factor two in 
our result, is only valid for our chosen spin-symmetric 
geometry. The investigation of more complex magnetic 
and structural geometries which might have lower con- 
ductance values will be the subject of future work. The 
role of fluctuations should also be considered, as they are 
likely to be important in a nanomagnetic system, partic- 
ularly at high temperature. 

IV. CONCLUSIONS 

We have generalized to magnetic transition metals 
the approach of Choi and Ihm for calculating the bal- 
listic conductance of an open quantum system. The 
method has been implemented with ultrasoft pseudopo- 
tentials and plane wave basis set in a DFT-LSDA ab- 
initio scheme. After checking the results against pre- 
vious calculations based on different methods, we ap- 
plied our method to calculate the ballistic conductance 
of monatomic Co and Ni nanowires with a single spin re- 
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versal. The s-like channel has a unit transmission prob- 
ability at all energies for both Co and Ni nanowires. On 
the contrary, the transmission probability of the heavy d 
electrons is heavily influenced by the spin reversal. For 
both wires we found that at the Fermi energy there is 
only one conductance channel per spin (s like) while all 
d electron ballistic conductance is completely blocked by 
the spin reversal, as w as p reviously guessed from band 
structure calculations |2g|. In the energy window be- 
tween 0.4 eV and 1.8 eV below the Fermi level the d 
electron conductance is smaller for the Co wire, due to 
an exchange splitting almost twice as large as that of 
the Ni wire. No structural distortions of linear chains 
(like zigzag structures or dimerizations) were considered 
in this paper, as they are not likely to be relevant in the 
short nanocontact connecting the two tips. More efforts 
to reproduce the realistic geometries of the nanocontacts 
will form the subject of future work. 
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